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Abstract. In phenomenological applications the time evolution of subsystems 
immersed in an external environment are sometimes described by Markovian 
semigroups of Redfield type that result non-positive: the appearence of negative 
probabilities is avoided by restricting the admissible initial conditions to those states 
that remain positive under the action of the dynamics. We show that this often adopted 
procedure may lead to physical inconsistencies in presence of entanglement. 
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1. Introduction 

The dissipative evolution of a system S immersed in a noisy environment E can be 
described via dynamical semigroups of linear maps j t , acting on density matrices 
p representing the states of the system. This reduced dynamics for S alone is 
obtained from the unitary time evolution of the full system S + E by tracing over the 
environment degrees of freedom and by further adopting a Markovian (i.e. memoryless) 
approximation. This procedure is physically justified when the interaction between S 
and E is weak and it has been successfully used in many phenomenological applications 
in quantum chemistry, quantum optics and atomic physics [1-12]. 

Nevertheless, the derivation of such reduced dynamics from the microscopic 
subsystem-environment interaction is often based on ad hoc, simplifying assumptions. 
As a consequence, the resulting reduced time evolutions are of Bio ch- Redfield type [7, 
8, 5] and therefore might not be fully consistent: typically, such naive dynamics do not 
preserve the positivity of the S density matrixjjl 

In order to cure these inconsistencies, the general attitude is to restrict the action of 
the non-positive dynamics to a subset of all possible initial density matrices, a procedure 
sometimes referred to as "slippage of the initial conditions" [13-16]. Physically, this 
prescription is ascribed to the short-time correlations in the environment, that are 
usually neglected in the derivation of the Markovian limit. 

In the following, we shall critically examine this prescription to cure inconsistencies 
produced by non-positive, Markovian evolutions and point out potential problems of 
this approach, in particular in presence of entanglement. More specifically, we shall 
study the behaviour of two subsystems, one immersed in the environment and evolving 
with a Redfield type dissipative dynamics, while the other does not evolve at all. We 
shall see that when the initial state of the two subsystems is entangled, the "slippage" 
prescription does not cure all possible inconsistencies of the two-subsystem dynamics. 

Preliminary, partial results on this line of investigation have been reported in 
Refs. [17-19]. In particular, in Ref. [18] the non-positive, dissipative evolution has been 
chosen in an ad hoc and rather abstract way in order to expose the just mentioned 
difficulties. Here, as reduced dynamics we adopt instead a Redfield non-positive 
evolution jt that has been used in various phenomenological applications [7-12]. With 
the help of both analytic and numerical methods, we shall then analyze the fate of the 
quantum correlations of two qubit systems when they evolve with a factorized dynamics 
7t <8) id: the "slippage" prescription is at work for the first qubit, while the other is inert 
and evolves with the identity operator. We shall explicitly show that in such a situation 
the purely local evolution 7$ ® id can increase the entanglement of the two systems, a 
clearly unphysical result. Therefore in presence of entanglement, the above mentioned 
prescription of restricting initial conditions to cure non-positive, Markovian dynamics 
does not seem completely satisfactory. 

| Exceptions to this general result are obtained using rigorous mathematical treatments [1-5]. 
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2. Single system dissipative dynamics 

We shall first study the dynamics of a single subsystem immersed in an external 
environment. As explained in the introductory remarks, the physical system we shall 
consider is a two-level system (qubit) immersed in a thermal bath. The system is 
described by 2 x 2 density matrices p, i.e. by positive Hermitian operators, with unit 
trace. On the other hand, the bath is modeled as an infinite dimensional reservoir in 
equilibrium at temperature T = f3~ l . Being infinitely large, the environment can be 
considered unaffected by the interaction with the subsystem and therefore to be in the 
reference equilibrium state 

-f3H E 

Tr(e~ 

where He is the hamiltonian describing the free dynamics of the bath. 

The total Hamiltonian describing the evolution of the compound system can be 
chosen as [8, 14-16] 

H = H s ®1 e + 1 s ®He + \Hi , (2) 

where 

H s = -a 3 , Hj = o x <g> B , (3) 

are respectively the subsystem Hamiltonian and interaction term, the parameter A 
playing the role of an adimensional coupling constant; at, i = 1,2,3 are the Pauli 
matrices and represent the subsystem operators, while B is an environment operator, 
taken for simplicity to satisfy the condition Tt e {peB) = 0. 

Using standard second order approximation in the coupling constant A and a naive 
Markovian limit, one finds that the time evolution of the reduced density matrix p for 
the system S is generated by a master equation of Bloch-Redfield type that takes the 
explicit form [7-9, 5]: 

^± = -i[H s ,p(t)} + L[p(t)} . (4) 

Besides the standard hamiltonian piece, the r.h.s. contains the extra contribution L, 
a linear map representing the effects of noise induced by the presence of the external 
bath. By setting, 

B(t) = e itHE B e~ itHE , (5) 
and further introducing the environment two-point correlation functions, 

G(t) = Tr[p E B(t)B] = Tr[p E BB(-t)] , (6) 
it can explicitly be written as [5]: 

L[p{t)} = \ 2 J ds^G{s) cos(o;s)[(Ti, <7ip(t)] — sin(o;s)[(T2, (Jip(t)] 

+ G(-s)\cos(ojs)[p{t)ai,Gi\ - sin(cjs)[p(t)o-i,o- 2 ]j [ . (7) 
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Therefore, the effects of dissipation and noise can be conveniently parametrized in terms 
of the following three phenomenological constants: 

roc 

a = X 2 ds cos(us)[G(s) + G(—s)] , 
Jo 

b = \ 2 ds sin(cjs) [G(s) + G(-s)] , (8) 

J 

d = i\ 2 ds sin(ws) \G(s) - G(-s)] . 
Jo 

Note that these parameters are not completely arbitrary. Indeed, since pE is a thermal 
state, it obeys the so called Kubo-Martin-Schwinger condition [20] : 

G(t) = G(-t - i/3) . (9) 

As a consequence, the parameters a and d above obey the following relation: 

a- d = e -/Jw (a + d) . (10) 

Further, one can show that the coefficient a must be positive [5]. 
The time evolution of the entries of the density matrix p, 

can now be explicitly given in terms of the parameters u, a, b and d: 
Pi(t) = l(l-^)(l-e- 2at )+pme- 2at , 

p 2 (t) = l- Pl (t) , (12) 
p 3 (t) = e- at {(cos(fit) -zl^±^sin(O^p 3 (0) + sin(fit) p*(0)} , 

where ft= [w 2 + 26cu - a 2 } 1 / 2 . 

Unfortunately, this evolution does not preserve the positivity of the eigenvalues of 
p for all times. In order to show this, it is sufficient to consider the following initial state 
p with entries: 



1A d\ A 1/ ,6\ 4a 2 -d 2 

* = 2V-to)> P ^A 1 + l a)i^W (13) 
It is a pure state, since Det[p] = 0. To have a p with positive spectrum, its determinant 
must remain non-negative for all times; in particular, its time derivative at t — must 
be positive, otherwise Det[p] would assume negative values as soon as t > 0. On the 
other hand, using (fl3l) . one easily sees that: 

d _ , „ N , a(46 2 + d 2 ) 

* Det <")i-=-i^W <0 - < 14 » 

a zero value being allowed only when b = d = 0. Further, note that because of the KMS 
condition (TlOl) . even if 6 = 0, a vanishing d can be obtained only at infinite temperature, 
i.e. when /3 = 0. It thus follows that at finite temperature the Markov approximation 
leading to the master equation (jlj), with generator as in ©, does not preserve positivity, 
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since a state like p is immediately turned into a matrix with negative eigenvalues as soon 
as t > 0. 

Although unphysical, the time evolution ( |T2l) generated by (BJ is nevertheless used 
in phenomenological applications because of its good asymptotic behaviour. In fact, it 
possesses a unique equilibrium state p eq , that can be easily determined by setting to 
zero the r.h.s. of (jl]): 



Therefore, it turns out that for asymptotically long times the system is driven to a 
thermal state at the bath temperature, a behaviour which is physically expected in such 
open quantum systems. 

In order to adopt the Redfield dynamics given in ffT2|) as a bona fide time evolution, 
a cure to the non-positivity needs to be introduced. The general solution that has 
been proposed is to restrict the space of initial conditions to those states p(0) that 
remain positive under the action of the Redfield dynamics. The general argument 
supporting this choice is that any Markovian approximation neglects a certain initial 
span of time, the transient, during which memory effects can not be ignored. During this 
short transient time, the environment acts in a very complicated way on the subsystem 
and the net result is the elimination of all states, like p in (Tl3|) . that would give rise to 
inconsistencies during the subsequent Markovian regime. This mechanism is known 
in the literature as "slippage of initial conditions" [13-16]. As we shall see in the 
next section, this prescription may cure the positivity preserving problem for a single 
subsystem, but appears to be inconclusive when dealing with bi- or multi-partite open 
systems in view of the existence of entangled states. 

3. Two qubit dynamics and entanglement 

We shall now extend the treatment discussed so far to the case of two qubits, one of 
which is still immersed in a heat bath and therefore evolves with the dissipative dynamics 
7t generated by the Redfield equations (jl]), while the other remains inert (it is usually 
called an ancilla) . The total time evolution for the two qubits is then in factorized form, 
7t <S> id, where "id" is the identity operator acting on the second qubit. In order to have 
a consistent time evolution, we shall further assume the "slippage prescription" at work 
for the first qubit: we remark that this prescription originates in the action of the bath 
during the transient and therefore can only involve the qubit inside the bath, and not 
the ancilla. 

Within this framework, we shall explicitly show the existence of states for the two 
qubits that 1) when traced over the ancilla degrees of freedom, belong to the set of 
admissible initial states for the non-positive dynamics j t , 2) remain positive under the 
action of the extended dynamics jt ® id and 3) nevertheless present an increase of their 
entanglement. This is clearly an unphysical result, because the evolution map acts 
locally, i.e. in a separate form, and therefore can not create quantum correlations. 




(15) 
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The existence of such states implies that the "slippage prescription" should take care 
not only of single system states developing negative eigenvalues, but also of possible 
inconsistencies related to the entanglement of these systems with any other ancilla. 

In order to explicitly expose this inconsistency, it will be sufficient to work within 
a special class of two-qubit density matrices, those for which the non- vanishing entries 
lie along the two diagonals: 

Pu\ 




P 





P22 
P23 





P23 
P33 



(16) 



/ Pn 


P23 P33 
\pl4 p 44 / 

Further restrictions on the entries of this matrix need to be imposed in order to represent 
a state. In particular, the trace must be one, P11+P22 + P33+P44 = 1, while the positivity 
of the spectrum implies the positivity of the two subdeterminants pupu — \pu\ 2 an d 
P22P33 _ I P23 1 2 and of the entries along the diagonal. The form fflBT) is particularly suited 
for our considerations since it is preserved by the action of the dynamics j t ® id; further, 
its entanglement content can be explicitly calculated. 

In this respect, a convenient measure of entanglement is provided by concurrence: 

C(p) = max{0, R 1 -R 2 -R 3 - Ra} , (17) 

where Ri are the square roots of the eigenvalues of R = p (02 <8> cr 2 ) p* (<7 2 <E> cr 2 ) taken 
in decreasing order; it vanishes for a separable state while takes positive values between 
zero and one for entangled states [21, 22]. For the state (jIB"|) . one explicitly finds: 

C(p) = max{0, max{|p 23 | - V 'P11P44, \pu\ ~ VP22P33}} • (18) 

It is then clear that the state ( fl6l) is entangled provided max{|p 2 3| — ^pnPu, 
\pu\ ~ v /p22P33} > 0. For simplicity, in the following we shall assume to start at t — 
with an entangled state p(0) fulfilling the more restrictive condition |p 2 3| — ^/pTIpIZ > 

I/flu I - VP22P33 > 0. 

Let us then consider the following two-parameter family of states 

I p 

- 3p) + |(1 - 2p) iv 

' -iv |^(3p-l) + |(l-2p) 





V 



-v 










p 



(19) 



where p and v are real constants satisfying the three constraints (necessary for positivity) 



1 -tf 



< p < 



l + i? 



3-2^ ^ 3 + 2^ 

-2 + 3?9 2 - V4 - 3^ 2 
9^ 



< p < 



-2 + 3d 2 + V4 - 3$ 2 
9^ 



(20) 



)v ;i_ 2 p)2-^2 (3/i _ 1) 2 >z/>;U> _ z/) 
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with $ = d/a. In writing (1201) . we have assumed a < b; this is not really restrictive since 
for a > b a similar family of states can be found. One can check that these matrices 
represent initially entangled two-qubit states that remain positive under the evolution 
7i (g> id; therefore, they are admissible states within the "slippage prescription" J§| 

However, these states present another, more subtle inconsistency than non- 
positivity. Indeed, using numerical methods, one can show that their concurrence 
increases for small times. The picture below displays the behaviour of the concurrence 
of one of these states as a function of timeljj] It shows an oscillatory behaviour that is in 
clear contradiction with quantum mechanics, since the dynamics is in factorized form. 



C[p] 




Figure 1. Concurrence behavior in time 



This unphysical behaviour of the concurrence can be studied analytically in the case 
of zero temperature; in fact, recalling ( flOl) . = implies the simplifying condition 
d = a. In this case, the matrix expression of our states and the corresponding constraints 
on the parameters p and v reduce to: 





/' 








-% v \ 







1-f/i 


iv 










—iv 


n 

2 





V 


-> 








n I 



(21) 



and 

0<M<jj, ~^2p-5p 2 >v> p> ^v . (22) 

With these conditions, it is easy to verify that p > 0. 

The evolution in time of this matrix under j t (g) id can be obtained from (Tl2j) ; using 
the labelling introduced in fTTol) for the entries of p, one explicitly finds: 

§ The system of inequalities (|20[) has solutions only when $ takes values in a certain range, which, 
recalling the condition (fT0|) , is related to the temperature of the bath. In our discussion we have taken 
V3/2 < <& < 1, since this allows certain simplifications in the calculations. 

|| The graph is drawn for the following representative values of the basic parameters: a/u) = 0.007, 
b/ui = 0.01, d/uj = 0.0065, with time in units of 1/uj. 
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Pu(t) 

P22{t) 

Pu(t) 



e~ 2at p , 

2 P 

P -2at 

2 

3p _ P e ~2at 

2 2 



Pu(t) = pn(t) 



P23(t) = f)23(t) 



-at 



-at 



av ,^ s bv . .av . . . 

-— cos(fit) - — sm(fit) +«— sin(f2t)(a; + 26) 
o ii o£2 



fi 



sin(fit) 



a 

T 



a; 



ii/ ( cos(fit) + — sin(fit) 



(23) 
(24) 

(25) 
(26) 
(27) 

(28) 



As previously mentioned, the positivity of the state p(t) at time t is assured by the 
positivity of the two sub-determinants pu(t)pu(t) — |pi4(t)| 2 and P22if) Pzzif) — I P23 ("^) | 2 ; 
in this case, these conditions read 

2 



2 v 



-2at\ 



p 3 

— 1 a 

2 V 2 P 



pe 



v 



-2at\ 



a _ . sin (fit) 
- cos(Ot) + & ^ ; 



+ a 2 



,sin (fit) / a; 



fi 2 V 



2 + 



> 



sin 2 (fit) fa 2 \ 2 / sin(fit)' 

' ( — + u + b) + cos(fit) + a- 



fi 2 V b 



fi 



(29) 



> . (30) 



In order to verify that these inequalities are indeed satisfied, recall from (jSJ) that a, b, 
d are proportional to A 2 ; since A is by assumption small, one can take a,l) <C w; being 



also fi 2 = uj 2 + 2uob 



u) 2 , we can neglect a and b with respect to uj and fi and 



then discard the terms proportional to a/fi, 6/fi and their powers with respect to those 
proportional toa/&ora>/fi~l. As a consequence, the conditions (29), (30) reduce to 

2a 2 



p 2 (3-e- 2at )>^-v 2 



P\ 1 - I^P-PZ 



-2at 



> 1v l 



(31) 

(32) 



these are easily seen to be satisfied thanks to the constraints in ( 1221) . In conclusion, 
the density matrices in fl2TT) are admissible initial states for the non-positive evolution 
7t ® id, since they remain positive for all times. 

Let us now compute their concurrence; one explicitly finds: 



C(p(t)) 



ue 



-at 



+ oo + b 



sin 2 (fit) 
fi 2 



cos(fit) + — sin(fit) 
P 



(33) 



- ^e- at VQ - 2e- 2at 
2 



(34) 



It is sufficient to examine the behaviour of C for small times: 
C(p(t))^u-p+^t + 0(t 2 ). 

Since a is positive, from this expression one immediately concludes that indeed C(p(t)) 
increases in time. 
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4. Discussion 

It is widely believed that the dynamics of a subsystem in weak interaction with an 
external environment can be described in terms of semigroups of linear maps 7 t generated 
by a Markovian master equation. In order to be physically acceptable, this effective 
description needs to satisfy basic physical requirements. In the first place, it must 
preserve the positivity of any initial density matrix, since their eigenvelues represent 
probabilities. In the second place, one has also to care of possible couplings with another 
system, not subjected to noise and inert, and therefore to guarantee the positivity- 
preserving character also of the semigroup of maps of the form j t ® id, as studied in the 
previous section. 

Unfortunately, most phenomenological derivation of reduced dissipative dynamics 
lead to semigroups of linear transformation that are not positive. To avoid 
inconsistencies, one usually restricts the possible initial states to those for which j t 
remains positive (the so-called "slippage of initial conditions"). This prescription works 
also in the case of the evolution ^ t <g> id for two subsystems, provided the initial state 
is in separable form: p(0) = J2iPiPi^ <8> pf \ Pi > 0, J2iPi = 1, where and pf 1 are 
admissible states for the first and second subsystems, respectively. 

However, as shown in the previous section, when the initial state p(0) is not in 
factorized form but still remains positive under the action of the non-positive dynamics 
7i <S> id, further, more subtle inconsistencies may arise. Indeed, we have found explicit 
examples of two-qubit states which under the action of j t ® id present an increasing 
concurrence. The creation of entanglement by a local operation is clearly unacceptable 
on physical grounds. This means that the "slipped" dynamics is still not free from 
inconsistencies. 

As a consequence, in order to continue to use non-positive reduced dynamics of 
Redfield type, a new, more general "slippage" mechanism should be invoked: it must 
take care not only of states developing negative eigenvalues but also of those presenting 
unphysical increase in entanglement. The only way to practically implement it is by 
further restricting the space of initial admissible states, discarding also some entangled 
ones. 

It should be noticed that possible inconsistencies are not limited to the two-qubits 
case; by considering more complicated ancillary coupling similar problems may arise for 
multipartite entangled states that should therefore also be eliminated by the "slippage 
operation". The risk of such a mechanism is to restrict too much the space of states, 
losing, in particular, many entangled states. These considerations seems to suggest 
that there is an intrinsic incompatibility between the existence of entangled states and 
the slippage prescription adopted to cure the inconsistencies that non-positive, reduced 
dynamics might produce. 

In closing, let us mention that in the few cases for which the Markovian limit of the 
subdynamics can be obtained in a rigorous way, the resulting evolution map j t turns out 
to be not only positive, but also completely positive [1-6]. In these cases, the compound 
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map 7 t <g> id is also completely positive and therefore no inconsistencies can arise, even 
when acting on entangled states. 
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